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Abstract 

Observations of the globular clusters NGC 6388 and M 15 were carried out by the 
H.E.S.S. array of Chcrenkov telescopes for a live time of 27.2 and 15.2 hours respectively. 
No gamma-ray signal is found at the nominal target position of NGC 6388 and M 15. 
In the primordial formation scenario, globular clusters are formed in a dark matter halo 
and dark matter could still be present in the baryon-dominated environment of globular 
clusters. This opens the possibility of observing a dark matter self-annihilation signal. 
The dark matter content of the globular clusters NGC 6388 and M 15 is modelled 
taking into account the astrophysical processes that can be expected to influence the 
dark matter distribution during the evolution of the globular cluster: the adiabatic 
contraction of dark matter by baryons, the adiabatic growth of a black hole in the 
dark matter halo and the kinetic heating of dark matter by stars. 95% confidence 
level exclusion limits on the dark matter particle velocity-weighted annihilation cross 
section are derived for these dark matter haloes. In the TeV range, the limits on the 
velocity-weighted annihilation cross section are derived at the 10~^^ cm^s~^ level and a 
few 10"^^ cm^s"^ for NGC 6388 and M 15 respectively. 

Subject headings: Gamma-rays : observations, Globular clusters. Black Holes, Dark 
Matter 
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Introduction 



Dark matter (DM) is expected to play a key role in the dynamics of a, wide range of sys- 
tems, from galactic object scale to galaxy cluster scale (jBertone et al. 20051 ). The inner parts 
of galaxies wer e used extensively to search for DM even though they are not necessarily dom- 
inated by DM (jEnglmaier &: Gerhard 20051 ). Observations by the H.E.S.S. experiment towards 
the Galactic Center reveal a very high energy ( VHE, > 100 GeV) g amma-ray signal with 



the flux level $(> ITeV) ~ 2 x lO'^^ (,j^-2g-i (jAharonian et al. 2009al ). The plausible DM- 



originated fraction of thi s signal is entirely dominated by standard astrophysical emission processes 
(jAharonian et al. 2006al ). As opposed to this, dwarf galaxies are believed to be am ong the most 
DM-dominated objects and are believed to have a reduced astrophysical background (jMateo 19981 ) 
since these objects have little or no recent star formation activity. Nearby dwarf galaxies in the 



Local Group ha ve already been observed with H.E.S.S 



I.e. 



tor and Carina (jAharonian et al. 20081 : jAharonian et al. 2009bl : jAharonian et al. 2011 



Sagittarius, Canis Ma 



or, Sculp- 
), yielding 

no signal. Exclusion limits on the velocity- weighed annihilation cross section between ~10~^^ to 



~10"22 cm^s^^ 



have been reported in the TeV range. 



Several Galactic Globular Clusters (GCs) have been observed with ground-based Cherenkov 
telescopes and upper limits on 7-ray emission from standar d astrophysical processes have been re- 
ported on Omega Ceiitauri, 47 Tucana, M 13, M 15 and M 5 (jKabuki et al. 20071 : Uharonian et al. 2009bl : 
Anderhub et al. 20091 : iMcCutcheon 20091 ) . GCs are also potential targets for indirect DM searches 
(jWood et al.|j2008j ). They are dense stellar systems of > 10 Gyr old, found in haloes of galaxies, with 
typical masses between 10^ and a few 10^ Mq, similar to dwarf galaxies. However, GCs are much 
more compact than dwarf galaxies. Observations of GCs do not suggest the presence of a significan t 
amount of DM, but rather that these objects a re dominated by baryons (jBinney &: Tremaine 19871 ). 
In the primordial formation scen ario of GCs (IPeebles 19841 1. GCs were formed in DM minihaloes 
before or during the reioni zation (IKomatsu et al. 20091 '). before formation of galaxies. However the 
distribution of GC colors ( Brodie h, Strader 20od ) suggests that only metal poor clusters have a 
cosmological origin while metal rich clusters formed in star-forming events such as galaxy-galaxy 
mergers. The existence of an ex tended dark halo require d by the GC primordial fo rmation scenario, 
has been challenged recently bv lBaumgardt et al.l (j2009l ) and lConrov et al.l (j2010l ). They show that 
the stellar kinematics of NGC 2419, a remote GC which experiences little tidal effe cts from the 



Milky Way, is incompatible with the presence of an extended dark halo. However, ICohen et al 



(j2010l ) have shown that the measured spread in the Ca abundance can be only explained if NGC 
2419 is the remnant of a more massive object. 

In the purpose of the paper, GCs are assumed to have formed in DM minihaloes and thus 
were DM-dominated in their primordial stage. Note that M 15 is a metal- p oor GC, [Fe/H]~- 
2.37 dHarris 199d ). while NGC 6388 is metal-rich, [Fe/H]~-0.55 (|Harris 1996l l. so the DM mini- 
halo scenario is better motivated for M 1 5 than NGC 6388. However, NGC 6388 might host a 
> 10^ MfD black hole (jLanzoni et al.l j2007j'l. Such massive (> 10^ Mq) black holes are not eas- 



ily formed in star-forming events (jYungelson et al. 20081 ) suggesting a primordial formation ori 
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gin. During the evolution of the globular cluster, the DM reacts to the infall of baryons and 
is pulled in towards the center. This process is usually referred to as the adiabatic contractioii 
(AC) model dBluinenthal et al. 1989 : IZeldovich et al. IQSol : Ijesseit et al. 2OO2I : ICnedin et al. 20041 : 
Prada et al. 2004 ). The effect of the contraction of DM in response to the baryon infall is partic- 
ularly important for th e calculation of the DM annihilation in bary onic environments such as the 
Galactic Center region (jCnedin &: Primack 2004 : IPrada et al. 20041 ) . 



The distr ibution of baryons and DM is affected by the kinetic heating o f DM by baryons 
(jMerrittl |2004| ) and by the presence of a black hole (BH) (jCondolo &: Silk 19991 ). A growing body 
of observations on GCs shows that they may harbor intermediate mass black holes (IMBHs) with 
masses ranging from 10^ to 10^ Mq, although the existence of these objects is not yet estab- 
lished. Indirect evidence includes the extrapolation of the empirical relation Meii-MBuige found 
for the sup ermassive BHs in gala ctic nuclei, which leads naturally to the prediction of existence 
of IMBHs ([Magorrian et al. 19981 ). Besides, ultra-luminous X-ray sources ( ULXs) apparently not 
associated to active galactic nuclei are p roposed to be related to IMBHs (IColbert &: Ptak 20021 : 
Swartz et al. 20041 : iDewangan et al. 20051 ). From the theoretical side, the existence of IMBHs is a 
generic assumption of scenarios that seek to explain the formation of supermassive BHs from the ac- 
cretion and merging of massive seeds (llslama et al. 20031: IVolonteri et al. 20031 : iKoushiappas et al. 20041). 
Several globular clusters may host IMB Hs: NGC 6388 (ILanzoni et al.ll2007l ) , u Centauri (jNovola et al. 20081 ) 
in the Milky Way or even Gl in M 31 dKong 200?! : luivestad et al. 200?! ). 



The paper is structured as follows. Section 2 describes the observations carried out with the 
H.E.S.S. experiment and the analysis of the data on the two Galactic GCs NGC 6388 and M 15. 
In section 3, the DM halo modelling relevant for the two GCs is briefly presented (more details 
are given in the appendix) and exclusion limits on the velocity-weighted annihilation cross section 
of DM particles are derived for realistic DM halo profiles. Finally, the results are summarized in 
section 4. 



2. Observations and data analysis 



The H.E.S.S. (High Energy Stereoscopic System) array of Cherenkov telescopes is located in 
the Khomas Highland of Namibia at an altitude of 1800 m above sea level. The system consists of 
four Imaging Atmospheric Cherenkov telescopes of 12 m diameter (~107m^) each. The total field 
of view of H.E.S.S. is 5° in diameter. The H.E.S.S. instrument achieves an angular resolution of 5' 
per gamma-ray and a point-like source sensitivity at the level o f 10~^^ cm~^s~^ above 1 Te V for a 



5(T detection in 25 hours at an observation zenith angle of 20° (|de Naurois &: Rolland 20091 ) 



For both M 15 and NGC 6388, data are taken in wobble mode (lAharonian et al. 2006bl ). where 
the pointing direction is ch osen at an alternating of fset of ±0.7° to ±1.0° from the target position. 
Standard quality selection (jAharonian et al. 2006bl ) is applied to the data. After the calibration of 
the raw data from PMT signals, the events are reconstructed using a technique based on a semi- 



- 6 - 



analytical shower development model (|de Naurois &: RoUand 20091 ). This reconstruction method 
yields a relative energy resolution of ~10% and an angular resolution of 0.06° per gamma-ray 
(68% containment rad i us). The backgrou nd level is estimated using the ring-background method 
([Piihlhofer et al. 20031 : iBerge et al. 20071 ). The core of the GCs is very small in comparison to 
the H.E.S.S. point spread function. They would thus be seen as point sources by H.E.S.S.. The 
source region referred to hereafter as the ON region, is defined as a circular region of 0.07°radius 
around the target position. The background region is defined as an annulus around the pointing 
position with a medium radius equal to the distance of the target and a half-thickness of 0.07°. It 
is referred to hereafter as the OFF region. The cut values used to select the gamma-ray events are 
shown in Tab. [TJ Given the Galactic latitude of the two Galactic GCs, contamination by diffuse 
TeV gamma-ray emission is very unlikely. The resu l ts presented below hav e been cross-checked 
by an independent analysis chain (|Berge et al. 20071 : lAharonian et al. 20051 ). Both analyses give 
consistent results. 



2.1. NGC 6388 



NGC 6388 is one of the best known Galactic GC. It is located at ~ 11.5 kpc from the Sun, at RA 
= 17^36^17.05' and Dec = -44°44'05.8" (J2000), and has a mass estimated to be ~ 10^ M©. The 
stellar mass density in the cor e, with radius = .4 pc (0.12 arcminutes), reaches '^5x10^ Mqpc"^. 
The tidal radius is rt ~ 25 pc (jLanzoni et al.ll2007l ). The structural properties of NGC 6388 are sum- 
marized in Tab.[2l Using the high-resolution HST and WFI observations at ESO. lLanzoni et al.l(|2007l ) 
show that the surface brightness density of stars significantly deviates from a flat core in the inner 
part, which is compatible with the existence of an IMBH with a mass of ~ 5 x 10^ Mq. A power law 
with a slope of -0.2 i s detected in the surface brightness density profile, which suggests the pr esence 
of a central IMBH teaumgardt et al. 20051 : iNovola fc Gebhardt 200fil : lUmbreit et al. 20081 ). The 
Chandra satellite has detected three X-ray sources, coincident in position with the centre of gravity 
of NGC 6388 locate d with an uncerta inty of 0.3". One of these may be the X-ray counterpart of 
the putative IMBH (|Nucita et al]2007l v 



The H.E.S.S. observations of NGC 6388 were taken between June 2008 and July 2009. The 
observation zenith angles range from 20° to 44° with a mean zenith angle of 22.9°, and the exposure 
time is 27.2 hours. The left-hand side of Fig. [T] shows the 6^ radial distribution of the ON and 
OFF events for NGC 6388. Nqn = 71 gamma-ray candidates were found in the ON region while 
the measured number of OFF events is Nqff = 1754. With the ratio of the off-source solid angle 
region to the on-source solid angle region a = 26.5, the expected number of background events 
in the ON region is 66.1. Since Nqn and Noff/^ are compatible at the 2a level, no significant 
gamma-ray excess is found above the background. An upper limit at 95% con fidence level (C.L .) 
on the number of gamma-rays can be derived using the method developed in iRolke et al.l(j2005l ): 

N95%C.L. ^ 21.6. 
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2.2. M 15 



M 15 (NGC 7078) is a well-studied Galactic GC centered at the position RA = 2li^28'^58.3'^ 
and Dec = 12°10'00.6" (J2000). It is situated at ~10kpc from the Sun. Its estimated mass is 
^ 5x 10^ M(7^. T he stellar mass density in the core with radius Tc = 0.04 pc is about 10^ Mqpc^'^ 
( Dull et alJl997l ). and the tidal radius is rt ~30pc. The structural properties of M 15 are summa- 
rized in Tab. [2j The surface brightness density of the GC M 15 sugg ests the presenc e of a stellar 



cusp in the inner part, at least down to d i stances of a few 10 ^ pc ( Dull et aljl997 ). M 15 



may 

thus harbor an IMBH ([Gerssen et al. 20021 : iKiselev et al. 20081 ) in its center. However, the study 
on n iilli-second pulsars in M 15 sets an upper limit of 10^ Mq on the mass of a hypothetical central 
BH (jDe Paolis et al. 19961 ). In what follows, no central back hole is assumed for the modelling of 
M 15. 

The observations of M 15 by H.E.S.S. were carried out in 2006 and 2007 with an offset angle 
of 0.7° and zenith angles from 34° to 44° resulting in 15.2 hours of high quality data at a mean 
zenith angle of 37.0°. The right-hand side of Fig. [T] shows the 0^ radial distribution of the ON and 
OFF events for M 15. The data analysis reveals no significant gamma-ray signal at the nominal 
position. With Nqn = 28, Nqff = 719 and a = 26.6, the expected number of background events 
in the ON region is 27.0. The 95% C.L. upper limit on the number of gamma-rays is _ 
11.5. Tab. [3] summarizes observations and upper limits on count numbers for NGC 6388 and M 15. 



3. Dark matter constraints 



3.1. Dark matter halo modelling 



The relaxation time, defined in Eq. (lAlj) of the appendix, has a muc h smaller value, Tr ^ lO'' yr , 
in GCs than in galaxies, where it is typically of the order of 10^^ yr ( Binney &: Tremaine 19871 ). 



^^^^^^^ 

core collapses ( 


Soitzer 


1987 


), the DM is comr 


(Zeldovich et al. 198nl: 


Bhimenthal et al. 19861 



This 

profil e is referred to hereafter as the AC NFW profile. But the kinetic heating of DM particles by 



stars (jMerrittll2004l ) tends to wash out the adiabatic contraction effect over a timescale of the order 



of Tf. Both effects were taken in to ac count in the modelling of M 15 and NGC 6388, following the 



approach of lMerritt et all (|2007l ) and lBertone k FairbarnI (j2008l ) 



As menti oned in the introduction, the primordial formation scenario of globular clusters 
(| Peebles 19841 ) assumed here requires that globular clusters were formed in extended DM haloes. 
The DM halo p rofile of a globular cl uster is thus modelled assuming an initial Navarro- Prenk- White 
(NFW) profile (jNavarro et al.lll997l ) described by: 



(1) 
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This DM halo is parameterized by a virial mas^ Mvir and a concentration parameter Cvir- The 
normahzation parameter pQ and the scale radius can be related to t he virial mass and the 
concentration parameter using the following relations (jNavarro et al.lll997l ) 



Cvir) Cvir 

where the function f{x) is, neglecting constants, the volume integral of the NFW profile given by 
/(x) = ln(l + x) — x/(l + x). The present baryonic mass of the GC provides a lower bound on 
its virial mass. Besides, for Mvir larger than 10^ M© the GC would be expected to spiral towards 
the center of the Milky Way in less than the ag e of the Universe . Conservative values for Mvir lie 
therefore in the range [5 X 10^ - 5 X 10^] Mq (|Wood et alfcoosi ). corresponding to Cvir between 
~48 and 65. In this paper, initial DM haloes of GCs are modelled with Mvir = lO'^ value 
of Cvir used in the model of NGC 63 88 is calculated fro m the formula of iBullock et al.l (j200ll ) . For 
M 15, the value of Cvir is taken from lWood et al.l(j2008l ). Both values are listed in Tab. HI 



The presence of a central BH changes the DM and stellar densities in regions where the 
BH dominates the gravitational potential, i.e. for distances to the BH lower than the radius of 
gravitational influence rjfl. The adiabatic growth of the BH leads to a spiked DM distribution with 
an index of 9/4 for an initial DM distribution with an index of 1, as for the NFW profile. This 
profile is referred to as the IMBH NFW profile. The spike is smoothed by the kinetic heating of 
DM by stars oyer th e timescale Tr, forming a density profile proportional to r~^/^ called DM crest 
(jMerritt et al.l 120071 ). which corresponds to the final profile. 



The dark halo models of M 15 and NGC 6388 ar e described i r i deta ils in the appendix. The 
DM halo of M 15 differs from the model published in lWood et al.l (|2008l ). since the effect of dark 
matter heating by stars is considered in addition to the effect of adiabatic contraction. The inferred 
DM mass densities of NGC 6388 and M 15 are shown in Fig. [2] and Fig. [3] respectively. 



3.2. Dark matter annihilation signal 



The gamma-ray fiux expected fr om DM annihilations can be decomposed into an astrophysical 
term and a particle physics term as ([Bertone et al. 20051 ): 

d^{An,E^) _ 1 {av) dN^ 
dK ~ 8^ ~^ 



Particle Physics 



X j{An)An . 

Astrophysics 



(3) 



^Mvir is defin ed as the mas s inside the radius Rvir assuming a mean density equal to 200 times the critical density 
of the Universe (lAmsler et al. 2 0091. 



^Mvir and Cvir are strongly correlated ( Navarro et al. 1997 . Bullock et al. 2001). In Bullock et al.l (I2OOI ) , Cy 



9 X (Mvir/1.5 X lO^^h ^Mq) where h is the present day normalized Hubble constant jAmsler et al. 20091 ). 
^The radius of gravitational influence of a BH is defined by the equation M(< rh) = JJ'' /5(r)d^r = 2Mbh 
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The astrophysical factor (J) is generally expressed as the integral over the line of sight (los) of the 
squared density averaged over the solid angle Ail: 

J=^l dn [ dsp\r{s)) (4) 
J An J los 

with r(s) = Y^s^ _|_ ^2 _ 2sso cos 9, sq the distance of the source from the Sun and 9 the opening 
angle of the integration cone centered on the target position. To match the analysis cuts used in 
this work (see Section 2), AQ is set to 5 x 10~^ sr. Tab. [4] shows the values of the astrophysical 
factor for the DM halo profiles presented in Sec. 13. li In the case of the IMBH NFW and final DM 
profiles for NGC 6388, the calculation of the astrophysical factor requires a minimum cutoff for 

—3/2 

the integration radius. For the IMBH NFW and final profiles, the integral diverges as r^^^ and 
log{r^l^) respectively, where rmm is the inner radius, rmm is usually taken as Max[r5,r^] where 
rs = 2GMbh / is the Schwarzschild radius of the black hole and ta is the self-annihilation radius 
calculated for an annihilation time of 10 Gyr. Typical values of niDM and (av) give rA — 10~^ pc 
so that rinin = ^A- The value of the astrophysical factor for the final profile is insensitive to the 
assumed value of r^i^. The values for M 15 are calculated using the DM modelling described in 
Sec. 13. 1[ The evolution of M 15 leads to a depletion of DM, implying a decrease of J. In the case 
of NGC 6388, the effect of the BH in the stellar environment boosts J to a value higher than that 
obtained for the initial NFW profile. 



3.3. Exclusion limits 



The exclusion limits on the particle physics parameters, i.e. the DM particle mass m-oM and 
the velocity- weighted annihilation cross section {crv), are calculated using the astrophysical factor 
calculated for each DM halo model by: 



{av) 



95% C.L. 
min 



87r 



m 



DM 



]y95% C.L. 



Tobs J (An) An 



(5) 



where dN^/dEy is the DM-produced diff^erential contin uum gamma-ray spectr um and ^eff is the 
effective area of the instrument during the observations (jAharonian et al. 20081 ). Fig. |4] shows the 
95% C.L. exclusion limits for NGC 6388 on (av) for the initial NFW (dashed line) and the final 
profile (solid line) deri ved from th e 95% C.L. upper limit on the number of gamma-rays. The generic 
parametrization from iBergstromI (|2000l ) as well as a para metrization including co ntributions from 
virtual internal Bremsstrahlung and fi nal state radiation (jBringmann et al. 20081 ) are used for the 
diff'erential gamma-ray spectra. The IBergstromI (j200d ) parametrization is derived from a fit to 
the gamma-ray spectrum from WIMP annihilations into W and Z pairs. The latter includes both 
gamm a-rays from virtual p articles and from charged particle final states of the pair annihilation of 
winos (jBertone et al. 20051 ). The limits are 1 to 3 orders of magnitude abov e the natural value o f 
the velocity- weighted annihilation cross section for thermally-produced DM (jBertone et al. 20051 ). 
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Fig. [5] shows the H.E.S.S. 95% C.L. exclusion limits for M 15 for the initial an d final DM profile s, 
as well as those obtained with the Whipple Cherenkov telescope (blue area) in IWood et al.l(|2008l ) . 
The thickness of the drawn lines represents the astrophysical uncertainty induced by the plausible 
mass range for the initial virial mass. The H.E.S.S. limit reaches {av) ~ 5 x 10~^^cm^s~^ and 
{av) 5 X 10"^'^ cm'^s"^ around moM = 2 TeV for the initial NFW profile and the final profile 
re spectively. For com parison, the exclusion limit obtained for H.E.S.S. using the DM halo modelling 
of IWood et al.l()2008l ) are also shown (gray area). Stronger constraints are obtained with H.E.S.S. 
due to the combination of larger effective area, 3 x 10^ m^ vs. 5 x 10^ m^, better angular resolution, 
0.06° vs. 0.25°, better upper limit on the number of gamma-rays (11.5 vs. ~ 67.7) and longer 
observation time (15.2 h vs. 1.2 h). 



4. Summary 

The present paper gives for the first time exclusion limits on dark matter towards several 
globular clusters taking into account all relevant astrophysical effects affecting the hypothetical 
DM halo. The H.E.S.S. observations reveal no significant gamma-ray excess from point-like sources 
located at the nominal position of the Galactic GCs NGC 6388 and M 15. The hypothetical DM 
halo has been modelled taking into account possible astrophysical processes leading to substantial 
changes in the initial DM profile: the adiabatic contraction of DM by baryons and the adiabatic 
growth of a BH at the center of the DM halo. The scattering of DM by stars in such a dense 
stellar environment has been taken into account to provide realistic final DM haloes. This effect is 
of crucial importance to model DM haloes in these baryon-dominated environments and leads to 
a depletion of DM during the evolution of the globular cluster. On the other hand, the presence 
of a central massive BH enhances the DM density in the center. The constraints on the velocity- 
weighted annihilation cross section of the DM particle are derived using DM halo profiles taking 
into account the above-mentioned astrophysical effects on the initial DM density. They lie at the 
level of a few 10~^^ cm'^s^^ for NGC 6388 in the TeV energy range. Assuming the absence of a 
massive BH in the center of M 15, the constraints are of the order of a few 10"^^ cm'^s"^. 
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Models of the M 15 and NGC 6388 dark matter halos 



A.l. M 15 dark matter halo 



The modelling of the M 15 DM halo proceeds in two steps. In the first step, the dark halo is 
assumed to be adiabatically compressed during the collapse o f the core of M 15. The model used for 
the initial and final baryon and DM densities is d escribed inlWood et al.l (|2008l ). The final baryon 
density is the observed mass density, taken from iGebhardt et al.l (|l997l ) . The DM is compressed 
during the ba ryonic collapse . The timescale for the collapse is ~ 100 Tj., where the relaxation time 
Tr is given by ISpitzeii (jl987l ) : 



3.4 X 10^ / Vj-ms 

In A Vkms^i 



m 



-2 



n 



pc 



-1 

-3 ' 



(Al) 



In Eq. (|Aip . Vrms is the velocity dispersion, n is the stellar density and In A is th e usual Couloni b 
logarithm. In the case of the center of M 15, taking Vrms = 10.2zfcl 4kms-^ (Dull et al]l997l ). 
In A = 13.1 and adopting a typical stellar mass value of m = 0.4 M0 (IDuU et al]l997l l. one finds 



Tr ~ 7x10^ yr. Tj. is an increasing function of the distance r to the center of the M 15. For 
f ^ '"heat = 5 pc, the relaxation time is larger than the age of the Universe. The central value of Tr 
and the position of r^eat have only weak dependencies on the actual values of Wrms and n when the 
latter are varied in their uncertainty ranges. Because of AC, the DM evolution takes place in less 
than a few orbital periods. The orbital period of a star orbiting the core of M 15 is of the order of 
1000 yr, which is much less than Tr. The AC method shoul d thus be valid. The val ue of the DM 



density after AC is similar to the result of lWood et al.l(l2008l ). see Fig. [3] and Fig. 7 of lWood et al. 
20081). At the same time, the DM is heated up by stellar matter. This process is described in 
MerrittI (j2004l ). DM is scattered by stars in a few Tr, leading to a depletion of the core. For r > 



5 pc, the DM distribution is not affected by heat i ng. T he DM scattering is taken into account with 
the procedure described in iBertone &: FairbarnI (j2008l ) for the M 4 globular cluster. A DM mass 
density of pMis ~ 35Mqpc^^ is obtained at the radius where the heating time is comparable to 
the age of the universe. For r < 5 pc, the DM halo is swept out by heating and the DM mass 
density was assumed to take the constant pMis value. The DM mass density of M 15 called final 
profile is shown on Fig. [3l 



A. 2. NGC 6388 dark matter halo 

The modelling of the NGC 6388 DM halo also proceeds in two steps. Th e first step is the AC 



of the DM halo by the IMBH and baryons. An initial baryon fraction of 20% (jSpergel et al. 20071 ) 
is assumed with the same spatial distribution like the DM. The AC scenario gives the resulting DM 
distribution knowing the measured baryonic mass profile. This DM halo profile is called AC NFW 
profile. The surface density profile of NGC 6388 is well fitted by a modified Ki ng model including a 
black hole, characterized by a core radius Vc = 7.2" and a concentration c = 1.8 ( Lanzoni et alj2007 ). 
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Using these parameters, the numerical integration of the Poisson equation yields the behavior of 
the gravitational potential from which it is straightforward to compute the baryonic density profile. 
The initial NFW profile is characterized by Mvir = 10^ M©. Since the mass of the IMBH is just 
a small fraction of the total mass in the core of NGC 6388, the dynamics of DM is influenced 
mainly by baryonic matter, except in the im mediate vicinity of the black hole. Using a central 
velocity dispersion of Urms = 18.9ib0.8 kms~^ ( Pryor &: Meylan I . ll993l ) and In A = 14.7, the central 
relaxation time is found to be ~ 8x10^ yr. The orbital period of a star orbiting the core of 
NGC 6388 is 5000 yr, so that the AC method is again valid. The relaxation time is larger than the 
age of the Universe for r > r^ea.t- The distribution of the DM density around the black hole (for 
r < rh) is changing with time, but tends to a power law with in dex 3/2 after a few T j- . The final 
DM distribution is thus obtained by extending the prescription of lBertone &: FairbarnI (j2008l ). Far 
from the center of the cluster, the stellar density is low and thus the heating time becomes large so 
that the DM distribution is unaffected. A mass density of ~140 Mqpc^^ is obtained at the radius 
''heat ~ 4 pc where the heating time is comparable to the age of the Universe. In the region between 
rh < r < rvipat,, the DM density is expected to be described by a smooth curve similar to Fig. 1 of 



Merritt et al.l ( l2007l ) . In the modelling for NGC 6388, the DM density was conservatively assumed 
to take a constant value of 14OM0pc~^ in the region < r < rheat- For r < rh, the DM density 
is given by ,o(r) = 14OM0pc"^ (r/rh)"^/^. The final profile for the DM distribution of NGC 6388 
is shown in Fig. [2j At the position of NGC 6388, the DM density from the smooth Galactic halo 
assuming a NFW profile is ~O.O3M0pc~^. 
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Table 1: List of H.E.S.S. analysis cuts ()de Naurois &: RoUand 20091 ). The shower depth is the 



reconstructed primary interaction depth of the particle and the nominal distance is the angular 
distance of the image barycenter to the camera center. 

Cut name 7-event cut value 

Shower goodness < 0.4 

Image charge (photo-electrons) > 120 
Reconstructed shower depth (radiation length) 

Reconstructed nominal distance (°) < 2 

Reconstructed event telescope multiplicity > 2 



This preprint was prepared with the AAS lATp^X macros v5.2. 
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Table 2: Properties of the NGC 6388 and M 15 globular clusters used in this study. 



NGC 6388 M 15 



Other Name / NGC 7078 

RA, Dec coordinates (J2000) 17''36™17.05^ -44°44'05.8" 2l'^28"'58.3", 12°10'00.6" 

Galactic coordinates (1, b) 345.55°, -6.73° 64.8°, -27.1° 

Distance (kpc) 11.5^ 10.0^ 

Core radius rc (pc) 0.4'' 0.04^ 

Tidal radius rt (pc) 25^ 30^= 

Estimated mass (Mq) 10^ 5x10^ 

Core stellar density (MqPc"^) 5x10^ lO'^ 

Central velocity dispersion Vrms (kms~^) 18.9±0.8^ 10.2±1.4'i 



Lanzoni et al.l(|2007h 



IPrvor fc Mevlan I (|l993 
iWood et al.1(|2008l ) 



Dull et al.l(|l997l ) 



Table 3: Exposure time, mean zenith observation angle, number of ON events, number of OFF 
events, a ratio and 95% C.L. upper limits on the number of gamma-rays from the H.E.S.S. obser- 
vations for NGC 6388 and M 15, respectively. See text for more details. 





NGC 6388 


M 15 


Exposure time^ (hours) 


27.2 


15.2 


Mean zenith observation angle 


22.9° 


37.0° 


NON 


75 


28 


NoFF 


1754 


719 


a 


26.5 


26.6 


]V^95%C.L. 


21.6 


11.5 



"After quality selection. 
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Tablc 4: Values of the los-integrated squared density averaged over the solid angle (J) expressed 
in units of 10^^ GeV^cm^^, for the different DM halo profiles. The integration solid angle is 
AQ = 5 X 10"*^ sr. The virial mass and concentration for the initial NFW profiles are given in 
brackets. 



DM halo profile name 


NGC 6388 


M 15 


Initial NFW (M^ir, Cyir) 


2.1 (lO'^ Mq, 60) 


1.5 (10"^ Mq, 50) 


AC NFW 


1.3x10^ 


4.3 xlO^ 


IMBH NFW 


2.2x10^ 


/ 


Final 


68 


14 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 



(a) NGC6388 (b) M15 

Fig. 1. — 9'^ radial distribution of the ON (gray histogram) and normalized OFF (black crosses) 
events for gamma-ray like events from a 0.07° radius region around the target position for (a) 
NGC 6388 (RA=17'^36™17.05" and Dec=-44°44'05.8", J2000) and (b) M 15 (RA = 2l'i28"^58.3^ 
and Dec = 12°10'00.6", J2000). The ON regions correspond to a maximum 6'^ of 0.005 deg^. No 
significant excess is found in the ON region for NGC 6388 or M 15. 
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Fig. 2. — DM and baryonic mass density distributions in NGC 6388. The DM density before (thick 
dashed hne) and after (thick dotted hne) the adiabatic contraction by baryons is shown. The initial 
DM distribution follows a NFW profile with Myir = lO'^'M©. The initial (thin dashed line) and 
final (thin dotted line) baryonic densities are displayed. The final DM density distribution after 
the effects of the adiabatic growth of the IMBH at the center of NGC 6388 and the kinetic heating 
by stars is presented (thick solid line). See text for more details. 
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Fig. 3. — DM and baryonic mass density distributions in M 15. The DM density before (thick 
dashed Hne) and after (thick dotted hne) the adiabatic contraction by baryons is shown. The 
initial DM distribution follows a NFW profile with Myir = IO'^Mq. The initial (thin dashed line) 
and final (thin dotted line) baryonic densities are displayed. The final DM density distribution 
after the effect of the kinetic heating by stars is presented (thick solid line). See text for more 
details. 
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Fig. 4. — H.E.S.S. upper limits at 95% C.L. on the velocity-weighted annihilation cross section (av) 
versus the DM mass ttibm for the Galactic GC NGC 6388. DM halo profiles shown here correspond 
to the initial NFW profile (dashed thick line) and the realistic profile taking into account plausible 
astrophysical effects (solid thick line). The contribution from internal Bremsstrahlung and final 
state radiation to the annihilation spectrum is also shown (dashed/solid thin lines) for both profiles. 
The natural value of (av) for thermally-produced DM is also displayed (long-dashed line). 
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Fig. 5. — H.E.S.S. upper limits at 95% C.L. on the velocity-weighted anihilation cross section {av) 
versus the DM mass ttidm for the Galactic GC M 15. Three DM haloes are shown: the initial NFW 
profile (dashe d thick/thin l ine), t he NFW profile after DM adiabatic contraction by baryons (AC 
NFW) used in I Wood et al. ( 2008 ) (gray area), and the final DM profile (solid thick/thin line). The 
effect of internal Bren isstrahlung i s also p resented (thin dashed/solid lines). The Whipple exclusion 
limits extracted from lWood et al.l ( 20081 ) is also plotted (blue area). The natural value of {av) for 
thermally-produced DM is shown (long-dashed line). 



